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Abstract 

We give evidence that the production of DCC's during a non- 
equihbrium phase transition can lead to an anomalous transverse dis- 
tribution of secondary pions when compared to a more conventional 
boost invariant hydrodynamic flow in local thermal equilibrium. Our 
results pertain to the linear a model, treated in leading order in large- 
A^, in a boost invariant approximation. We also show that the interpo- 
lating number density of the field theory calculation plays the role of 
a classical relativistic phase space number distribution in determining 
the momentum distribution of pions in the center of mass frame. 



1 Introduction 



Recently there has been a growing interest in the possibility of producing 
disoriented chiral condensates (DCC's) in a high energy collision |l|, |]. 
This idea was first proposed to explain CENTAURO events in cosmic ray 
emperiments where there was a deficit of neutral pions 0]. It was proposed 
that a nonequilibrium chiral phase transition such as a quench might lead to 
regions of DCC 0. 

In a previous paper 0, we discussed in detail the time evolution of a 
non-equilibrium chiral phase transition in the linear a model in the leading 
order in large-N approximation. We showed in that paper that the expansion 
into the vacuum of the initial energy distribution led to rapid cooling. This 
caused the system, initially in quasi local thermal equilibrium to progress 
from the unbroken chiral symmetry phase to the broken symmetry phase 
vacuum. This expansion is accompanied by the exponential growth of low 
momentum modes for short periods of proper time for a range of initial condi- 
tions. This exponential growth of long wave length modes is the mechanism 
for the production of disordered chiral condensates. Thus the production of 
DCC's results in an enhancement of particle production in the low momen- 
tum domain. Whether such an instability occurs depends on the size of the 
initial fluctuation from the initial thermal distribution. The relevant mo- 
menta for which this exponential growth occurs are the transverse momenta 
and the momenta kr, = —Ez + tp conjugate to the fluid rapidity variable 
7] = iaiihr^{z/t). We also found that the distribution of particles in these 
momenta had more length scales than found in local thermal equilibrium. 
When there is local thermal equilibrium, the length scales are the mass of 
the pion and the temperature which is related to the changing energy density, 
both of which depend on the proper time r. 

In this paper we reexpress our previous results in terms of the physically 
measurable transverse distribution of particles in the collision center of mass 
frame. We show that whenever DCC's are produced there is a noticable dis- 
tortion of the transverse spectrum, namely an enhancement of particles at 
low transverse momentum, when compared to a local equilibrium evolution. 
We will consider two cases, one in which there is exponential growth of low 
momentum modes due to the effective pion mass going negative during the 
expansion, and one where the initial fluctuations do not lead to this exponen- 
tial growth. Both situations will be compared to a purely hydrodynamical 
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boost invariant calculation based on local thermal equilibrium. The idea 
that the production of DCC's might lead to a distortion of the transverse 
distribution of pions is also found in the work of Gavin and collaborators 
[P]. However, our calculation is more realistic in that we exactly solve the 
field theory in a large approximation in the boost invariant approximation, 
using the expansion of the plasma as the cooling mechanism for producing in- 
stabilities. In determining the actual spectra of secondaries, we find that the 
adiabatic number operator of our large-N calculation replaces the relativistic 
phase space density g{x,p) of classical transport theory in determining the 
distribution of pariticles in rapidity and transverse momentum. This makes 
it easy to compare our results with a standard hydrodynamical calculation 
in the boost invariant approximation which assumes the final pions are in 
local thermal equilibrium in the comoving frame. 

The model we use to discuss the chiral phase transition is the linear sigma 
model described by the Lagrangian: 

L = ^d<l>-d<l>-^X{<^-<^-v^f + Ha. (1.1) 

The mesons form an 0(4) vector $ = (a, tTj) This can be written in an 
alternative form by introducing the composite field: x = " ^ ~ ^^)- 

U = + X)<l^^ + ^ + \xv' + Ha (1.2) 

The effective action to leading order in large is given by 

x]= j rf'x[L2($, X, H) + ^iVtr IuGq^] (1.3) 

G,\x,y) = ip + x{x)] 5\x-y) 
Varying the action we obtain: 

[□ + X(a:)]7r, = [□ + x(x)]a = (1.4) 

where here and in what follows, 7rj,cr and x refer to expectation values. Vary- 
ing the action we obtain 

X = -Xv^ + X{a^ + vr ■ tt) + XNGo{x, x). (1.5) 
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If we assume boost invariant kinematics which result in flat rapidity dis- 
tributions, then the expectation value of the energy density is only a function 
of the proper time. The natural coordinates for boost invariant {v = z/t) 
hydrodynamical flow are the fluid proper time r and the fluid rapidity 77 
defined as 

To implement boost invariance we assume that mean (expectation) values of 
the fields $ and x are functions of r only. We then get the equations: 

T-^dr rdr <f.(r) + x{r) <f.(r) = H6,^ 

^{r) = \{-v^ + <|.2(r) + iVGo(x, x; r, r)), (1.6) 

To determine the Green's function Go{x,y;T,T') we introduce the auxiliary 
quantum field 0(x, r) which obeys the sourceless equation: 



[r-'drTdr -r-^dl-dl + x{x))<p{x,T) = 0. (1.7) 

Go{x, y- r, r') =< T{</)(x, r) </)(y, r')} > . 
We expand the quantum fields in an orthonormal basis: 

<^(^?,a;_L,r) = j [d^k](exp(?fca;)/k(r) aj, + h.c.) 

where kx = k^r] + k±x±, [rf^k] = dkr^(Pk±/ (27r)^. The mode functions and x 
obey: 

/k + (^ + + X(r) + ^)/k = 0. (1.8) 

x{r) = \[-v' + e,{T) + ^N j[d\]\Ur)\'' {l + 2n^)). (1.9) 
We notice that when x goes negative, the low momentum modes with 

+ 1/4 



r2 



+ ki< \x\ 



grow exponentially. However these modes then feed back into the x equation 
and this exponential growth then gets damped. It is these growing modes 
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that lead to the possibhty of growing domains of DCC's as well as a modifica- 
tion of the low momentum distribution of particles from a thermal one. The 
parameters of the model are fixed by physical data. The PCAC condition is 

In the vacuum state Xo<^o = ^n'^o = H, so that ctq = 92.5 MeV. The 
vacuum gap equation is 



Jo 



'0 2^k^ + ml 
This leads to the mass renormalized gap equation: 
x(t) - ml = -\fi + AtKr) 

+Aa. I [d\]mr)? (1 + 2 ,,,) - ^^J^^J . (1.10) 

A is chosen to fit low energy scattering data as discussed in 0. 

If we assume that the inititially (at tq = 1) the system is in local thermal 
equilibrium in a comoving frame we have 



where /?o = I/Tq and = J^ + kl + ^(ro). 



'0 



The initial value of x is determined by the equilibrium gap equation for 
an initial temperature of 200MeV and is .Ifirr"^ and the initial value of a is 
just The phase transition in this model occurs at a critical temperature 
of lQ(]MeV . To get into the unstable domain, we then introduce fluctuations 
in the time derivative of the classical field. We varied the value of the initial 
proper time derivative of the sigma field expectation value and found that 
for To = 1/m there is a narrow range of initial values that lead to the growth 
of instabilities, namely .25 < \a\ < 1.3. 

Figure one displays the results of the numerical simulation for the evo- 



lution of the system (|1.8| )- ([L9| ). We display the auxiliary field x units 
of fmr"^ and the proper time in units of fm {lfm~^ = 197MeV) for two 
simulations, one with an instability {ctIto = ~1) and one without {crlro = 0). 
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Wc notice that for both initial conditions, the system eventually settles down 
to the broken symmetry vacuum result as a result of the expansion. 

To determine the spectrum of particles we introduce the interpolating 
number density which is defined by expanding the fields in terms of mode 
functions /° which are first order in an adiabatic expansion of the mode 
equation. 

= dvk/dt^Uk, (1.12) 

where a;k(T) = (/c^/t^ + A;i + x(T))^''^. This leads to the alternative expansion 
of the fields: 

0(77, r) = :^J m{eM^kx)f^{r) a^{r) + h.c.) (1.13) 

The two sets of creation and annihilation operators are connected by a 
Bogoliubov transformation: 

ttkir) = a{k,T)ak + P{k,T)a^_f.. (1-14) 
a and (3 can be determined from the exact time evolving mode functions via: 

o^iKr) = lUk -^--^/fe) 

Pik.r) = ^ifk^-^fk). (1.15) 
In terms of the initial distribution of particles no{k) and (3 we have: 

nk{r) = f{k^,k±,T) =< 4(r)afc(r) > 

= no{k) + \p{k,T)Wl + 2noik)). (1.16) 

nfc(r) is the adiabatic invariant interpolating phase space number density 
which becomes the actual particle number density when interactions have 
ceased. When this happens the distribution of particles is 

= ir^dxldkldr^dk, ^^-^^^ 

We now need to relate this quantity to the physical spectra of particles mea- 
sured in the lab. At late r our system relaxes to the vacuum and % becomes 
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the square of the physical pion mass m . We introduce the outgoing pion 



particle rapidity y and m± = \Jkj_ + m"^ defined by the particle 4-momentum 
in the center of mass coordinate system 

— {m±coshy, k±,'m±smhy) (1-18) 

The boost that takes one from the center of mass coordinates to the comoving 
frame where the energy momentum tensor is diagonal is given by tanhrj — 
V = z/t, so that one can define the "fiuid" 4- velocity in the center of mass 
frame as 

m'' = (cosh ?7, 0,0, sinhr/) (1-19) 
We then find that the variable 



(^k = \lml + § = (1.20) 



has the meaning of the energy of the particle in the comoving frame. The 
momenta A;^ that enters into the adiabatic phase space number density is 
one of two canoncial momenta to the variables defined by the coordinate 
transformation to fiuid light cone variables. Namely the variables 



2 \t-z 
have as their canonical momenta 

kr = Et/r - k^zjr A;^ = -Ez + tk^. (1.21) 

To show this we consider the metric ds^ — dT^—T^dif' and the free Lagrangian 

m dx^ dx^ 

= — 9ui/ (1-22) 

2 ^^'^ ds ds ^ ' 

Then we obtain for example 

dr r,dr, dt ,dT , dz. 



ds dt ds dz ds 
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The interpolating phase-space density / of particles depends on A;^, k_L, r, 



and it is found to be ?7-independent and given by eq. (|1.17|) . 

In order to obtain the physical particle rapidity and transverse momentum 
distribution, we change variables from (77, A;^) to {z, y) at a fixed r where y is 
the particle rapidity defined by (|1.18|) . We have 



E- 



J ndz dx\ J f{kn, r) (1-24) 



d^k ndydk^ 
where the Jacobian J is evaluated at a fixed proper time r 

J -- 



dkri/dy dkri/dz 




dkrf drj 


drj /dy drj/ dz 




dy dz 


m_L cosh(?7 — y) 


dkri 


1 T ' 


coshr] 


dz 





We also have 



m_L cosh(r7 — y); k^f = —Tm± sinh(r7 — y) . 



:i.25) 



:i.26) 



Calling the integration over the transverse dimension the effective transverse 
size of the colliding ions A± we then obtain that: 



d^N 
Trdy dk'j_ 



A_L J dknf{kr^,k±,T) 



[1.27) 



This quantity is independent of y which is a consequence of the assumed 
boost invariance. 

In a hydrodynamical model of heavy ion collisions 0, the final spectra of 
pions is given by a combination of the fluid flow and a local thermal equilib- 
rium distribution in the comoving frame. One calculates this spectra at the 
critical temperature Tc{x, t) when the energy density goes below ec = (Yi/mc)^ ■ 
This defines the breakup surface, after which the particles no longer interact 
so that this distibution is frozen at that temperature. For an ultrarelativistic 
gas of pions, this occurs when ~ m. The covariant form for the spectra of 
particles is given by the Cooper- Frye formula 



E 



d^N d^N 



d^k TTdkj_dy 



g{x, k)k^daf, 



[1.28) 
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Here g{x, k) is the single particle relativistic phase space distribution func- 
tion. When there is local thermal equilibrium of pions at a comoving tem- 
perature T'c(r) one has 

9{x, k) = ^{exp[P«^/Tj - 1}-^ (1.29) 

In the boost invariant approximation, the temperature is only a function of 
the proper time r and the break up surface is a constant r surface: 

da^ = A^{dz,0,0,dt) = A ^dr] {cosh r],0,0,smhri) (1.30) 

where Aj_ is the effective transverse size of the system following the heavy 
ion collision. Therefore we obtain the relationship that 

k^da^ = A j_mj_T cosher] — y) = A^\dkri\ (1-31) 

Thus we can rewrite our expression for the field theory particle spectra as 

d^N 



Tidy dk] 



A± J dkjikr,, k^, t)= J fikr,, k^, T)k>^da^ (1.32) 



where in the second integration we keep y and r fixed. This shows that the 
Cooper-Frye formula for the final spectra of secondaries, which was derived 
from classical transport theory, also applies in field theory when one has a 
gaussian density matrix as in the large-N approximation. It also reconfirms 
the idea that the interpolating phase space number density plays the role 
of a classical transport phase space density function, as was found in our 
calculation of pair production from strong electric fields ||^ 

In Figures 2 and 3 we compare the boost invariant hydro dynamical re- 
sult for the transverse momentum distribution using critical temperatures 
of Tc = 140, 200 MeV to the two nonequilibrium cases represented in figure 
1. Figure 2 pertains to the initial condition a\ro = —1 In this case there is 
a regime where the effective mass becomes negative and we see a noticable 
enhancement of the low transverse momentum spectra. We have normalized 
both results to give the same total center of mass energy Ecm- Figure 3 
corresponds to the initial condition a\ro = 0. Here we notice that there is a 
little enhancement at low transverse momenta. 
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In conclusion, we see that a non-equilibrium phase transition taking place 
during a time evolving quark-gluon or hadronic plasma can lead to an en- 
hancement of the low momentum transverse momentum distribution. In 
particular, if a Centauro type event is not accompanied by such an enhance- 
ment one would be suspicious of ascribing this event to the production of 
disoriented chiral condensates as a result of a rapid quench. 
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FIGURE 1. Proper time evolution of the x fi^ld for two different initial 
conditions ((t(1) = and ct(1) = -1) with = 92.5MeV". 

FIGURE 2. Single particle distribution k^^^ vs. k± in the mean field 
approximation for the initial condition = —1 compared to two different 
hydrodynamical calculations with two different breakup temperatures: Tc = 
rTT-Tr and Tc = lAmT^. We have normalized all the curves to give the same 
center of mass energy. 
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FIGURE 3. Same comparison as in fig. 2 but with the mean field ap- 
proximation for the initial condition cr(l) = 0. 
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da/dT=0 initial conditions 
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